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Abstract. In this paper the number and lengths of minimal length lattice knots 
confined to slabs of width L, is determined. Our data on minimal length verify 
the results by Ishihara et. al. [8] for the similar problem, expect in a single case, 
where an improvement is found. From our data we construct two models of grafted 
knotted ring polymers squeezed between hard walls, or by an external force. In 
each model, we determine the entropic forces arising when the lattice polygon is 
squeezed by externally applied forces. The profile of forces and compressibility of 
several knot types are presented and compared, and in addition, the total work 
done on the lattice knots when it is squeezed to a minimal state is determined. 
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1. Introduction 

Chemically identical ring polymers may be knotted and these are examples of 
topological isomers which may have chemical and physical properties determined by 
their topology. There has been a sustained interest in the effects of knotting and 
entanglement in polymer physics and chemistry, and it is known that entanglements 
may play an important role in the chemistry and biological function of DNA 29J . For 
example, entanglement and knotting are active aspects of the functioning of DNA and 
are mediated by topoisomerases [161 117) , while proteins are apparently rarely knotted 
in their natural active state [2"D] . 

Ring polymers with specified knot type have been chemically synthesized [S], 
but more often, random knotting of ring polymers occurs in ring closure reactions 
[23] [2j [20] . In this case, a spectrum of knot types are encountered EH EI] , and 
these are a function of the length of the polymer: Numerical studies show that longer 
ring polymers are knotted with higher frequency and complexity |14) . 

Ring polymers adsorbing in a plane or compressed in a slab also appear to have 
increasing knot probability |19[ 121) . although the probability may decrease in very 
narrow slabs |2 7) . Similar effects are seen when a force squeezes a ring polymer in a slab 
[ID] , and the results of the calculation in reference [DJ suggest that knotted polygons 
will exert higher entropic forces on the walls of a confining slit. More generally, the 
phase behaviour of lattice ring polymers confined to slabs and subjected to external 
forces were examined in reference [25] . 

The entropic force of a knotted ring polymer confined to a slab between two plates 
were examined using a bead-spring model in reference [18^ . In this study it was found 
that more complex knot-types in a ring polymer exert higher forces on the confining 
walls of the slab (if the slab is narrow). 

In this paper we obtain qualitative results on the entropic properties of tightly 
knotted polymers confined to a slab or squeezed by a flexible membrane, using minimal 
length cubic lattice knots. We will consider two different models. 

The first is a model of a tightly knotted ring polymer of fixed length squeezed 
between two hard walls or plates (see figure [T]). Self- avoidance introduces steric 
repulsions between monomers which causes (self)-entanglemcnt of the polymer. 
Confining the polymer to a slab results in the loss of configurational entropy inducing 
a repulsive force which depends on the entanglements between the walls of the slab. 
This entropically induced repulsion will be overcome by an externally applied force at 
critical magnitudes, and we shall determine these critical forces for several knot types 
in our model. An external force will tend to squeeze the walls of the slab together, 
and at some critical widths of the slab, the polymer cannot shrink further without 
expanding laterally and increasing its length. Beyond this critical width there may 
also be an elastic energy contribution to the free energy - the polymer stretches in 
length to accommodate the narrow slab. We model this with a Hooke energy. 

The second model is inspired by the study in reference [7]. In figure 1 therein, 
a polymer is grafted to a hard wall and covered by a soft flexible membrane. The 
membrane may be modeled by a hard wall as in the first model above. On the other 
hand, a pressure difference across the membrane will induce forces on the monomers 
in the top layer of the grafted polymer. This couples the force to the height of the 
monomers in the top layer - the force is mediated directed through the membrane and 
push on the highest monomers of the polymer. 

We consider the models presented above and illustrated in figures [I] and [2] in turn. 
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Figure 1. A minimal length lattice knot of type 3i squeezed between two hard 
walls. The lattice knot is grafted to the bottom wall - that is, it has at least one 
vertex in this wall. If the two hard walls are a distance L apart, then the height 
of the lattice knot is h < L, since the polygon must fit between the two walls. 
The number of these lattice knots of length n are denoted p^(K). 

In section 2 we present our models and discuss the collection of numerical data. An 
implementation of the GAS-algorithm for knotted cubic lattice polygons in slabs [TT] 
was used, and we collected data on the entropy and minimal length of knotted lattice 
polygons up to eight crossings. Our data verify similar results obtained in reference 
[5]. In section 3 we discuss the first model and present our data, and in section 4 we 
present data on the second model and discuss our results. The paper is concluded in 
section 5. 




Figure 2. A minimal length lattice knot of type 3i grafted to the bottom wall 
and covered by a weightless flexible membrane. The lattice knot is grafted to the 
bottom wall (which is hard) by having at least one vertex in this wall. A pressure 
gradient over the top membrane will exert a force compressing the knot onto the 
bottom wall. The force is mediated through the membrane onto the monomers 
in the top layer of the lattice knot. If the height of the lattice knot is h, then 
the number of vertices at height h is conjugate to a force pushing the lattice knot 
towards the bottom wall. The number of these lattice knots of length n is denoted 
by Pn(K; h). 
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2. Models of lattice knots in slabs 

A lattice polygon is a sequence aj n = {vo,vi, . . . ,v n -i} of distinct vertices Vi such 
that ViVi+i, for i = 0, 1, . . . , n — 2 and v n -ivo are unit length edges in the cubic lattice 
Z 3 . Two polygons loq and u>\ are equivalent if they are translates of each other. The 
length of a polygon is its number of edges (or steps). A polygon is, by inclusion, a 
tame embedding of the unit circle in M 3 and so has a well defined knot type. A lattice 
polygon with specified knot type is a lattice knot. 

The number of distinct polygons of length n and knot type K will be denoted 
p n (K). For example P4(0i) = 3 while p n (0i) = if n < 4 or if n is odd, where 
0i is the unknot in standard knot notation. It is also known that j>24(3^) = 1664 
while p n (3i ) = if n < 24 [22], where 3^ is the trefoil knot type. By symmetry, 
p n (K + ) — p n (K~) if AT is a chiral knot type. 

The minimum length of a lattice knot K is the minimum number of edges required 
to realise it as a polygon in Z 3 . For example, the minimal length of knot type 0i is 4 
and of knot type 3^ or 3|f is 24 edges. The minimal length is denoted by uk, so that 
n 0l = 4 and n 3 + — 24 4 . It also known that n 4l = 30 (4i is the figure eight knot) 
and n 5 + = 34 [22]. Beyond these, only upper bounds on tir are known. 

If v € Lu n is a vertex in a lattice polygon, then the Cartesian coordinates of v are 
(X(v),Y(v),Z(v)). A polygon is grafted to the XY plane Z = if Z(v) > for all 
vertices v in ui, and there exists one vertex, say Vq, such that Z(vq) = 0. The height 
of a grafted polygon is h = max {Z(i>) | v 6 

A grafted polygon ui n is said to be confined to a slab of width L if < Z(u) < i 
for each vertex v G w„. A polygon confined to a slab is illustrated in figure Q] where 
the bottom wall and top waH of the slab are indicated. 

Next we define the number of polygons of length n, knot type K, and confined 
in a slab of width L by p^(K). Similarly, define the number of polygons of length n, 
knot type K, and of height h by p n (K; h). Clearly, 



The minimal length of a lattice knot in a slab of width L will be denoted ul^k- 
For example, one may deduce that n2,3 1 = 24 from reference [4]. However, simulations 
show that = 26 [8]. In other words, in a slab of width L = 1, it is necessary to 
have a polygon of length 26 to tie a lattice trefoil, while if L > 2 then 24 edges will 
be sufficient (and necessary) [4]. 

With these definitions, we define two models of grafted lattice knots. The first 
model is that of a grafted lattice knot in a slab with hard walls (see figure [lj . In 
this model the confinement of the lattice knot will decrease its entropy, and this will 
induce an entropic repulsion between the top and bottom walls of the slab. The 
discrete geometry in this model implies that the induced entropic force is given by 
free energy differences if the distance between the walls is decreased by one step. 

The second model is illustrated in figure [2] The grafted lattice knot is covered by 
a flexible and weightless membrane. A (positive) pressure difference in the fluid above 
and below the membrane induces a force pushing on the top vertices in the lattice 
knot. A negative pressure difference in the fluid results in an effective pulling force on 
the vertices in the top layer of the lattice knot. In this model, the partition function is 
given by all the states of grafted lattice knot (including those of any height). The force 
/ induced by the pressure difference will push on the highest vertices in the polygon 




(1) 



h<L 
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as it is mediated by the flexible lightweight membrane onto these vertices. A (linear) 
compressibility of the lattice knot can be determined by taking second derivatives of 
the free energy of this model to the applied force, as we shall show below. 

2.1. Numerical approach 

In this paper we examine the properties of minimal length lattice knots squeezed 
between two hard walls, or squeezed by an applied force towards a hard wall. In both 
these models it is necessary to determine the number and length of lattice knots in 
slabs of width L. Some data of this kind were obtained in reference jS], and we will 
at the same time verify in most cases, and improve in one case, on their results. 

Our numerical approach will be the implementation of the GAS-algorithm for 
lattice knots [TTJ [T2J [T3] using BFACF elementary moves [TJ |3] ■ The lattice knots will 
be confined to slab S^ of width L defined by S^ =1? X {0, 1, . . . , L}. We implement 
the algorithm by noting that BFACF elementary moves on unrooted cubic lattice 
polygons are known to have irreducibility classes which are the knot types of the 
polygons. The proof of this fact can be found in reference [15]. Note that the proof in 
reference |15) applies mutatis mutandis to the model in this paper as well, provided 
that L > 2. 

We estimate p nL k (K; h) (the number of polygons of knot type K, height h and 
of length ul.k) using the GAS algorithm for knotted polygons. Here n^jf is the 
minimal length of grafted lattice knots of type if in a slab of width L. By summing 
h < L, one obtains (K), the number of grafted lattice knots of type K which can 
fit in a slab of width L, of length n,L,K (and thus of height h < L). Our results are not 
rigorous, and in a strict sense the results of til,k are upper bounds while p riL K (K; h) 
are lower bounds. Since the GAS algorithm can be implemented as a flat histogram 
method, it is efficient at rare event sampling and thus at finding knotted polygons of 
minimal length. A comparison of our results with the data in reference [8] makes us 
confident that our results are exact in almost all cases. 

Data on lattice knots in slabs Sj,, with L > 2 can be collected as in reference 
[12] . The case L = 1 requires further scrutiny. Data in this ensemble were collected 
by generating lattice knots in §2, and sieving out lattice polygons which fit into Si. 
By biasing the sampling to favour polygons which fit into Si, we were successful in 
generating lists of lattice knots in Si . We are reasonably confident that in most cases 
our lists of knotted polygons are complete. 

We display our data on the minimal length of lattice knots in Si, in appendix A 
in tables [5] and [7] (for some compound knots) . The data in table [5] agree for all knot 
types, except for 8is, with the data in reference [8]. We improved on the estimate for 
the minimal length of 818 in Si by finding states at length n = 70, compared to 72 in 
that reference. 

As one might expect, we observe a steady increase of rij,^ with decreasing L 
for each knot type, and also down table [SJ Decreasing L squeezes lattice knots in 
narrower slabs, and at critical values of L there is an increase in the minimal length. 
For example, for the trefoil knot 3^, there are realizations of lattice knots with this 
knot type at n = 24 edges for L > 2. However, if L = 1, then 26 edges are needed. In 
figure [3] the minimal lengths of lattice knots in Si are compared to the minimal lengths 
of lattice knots in the bulk lattice. In this scatter plot each knot type has coordinates 
(nx t R, tik)- We found that n\^x > na for all non-trivial knot types, but the data do 
cluster along a line showing a strong correlation between these two quantities. For 
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Figure 3. A scatter plot of minimal lengths n^jf in Si (X-axis) and minimal 
lengths nj{ (y-axis). 



slabs Sl our data are displayed in tabled showing that ri\ K > nL,K > nx generally. 

Naturally, til k is a non- increasing function of L for a given knot type K . In some 
cases there is a large increase in til k with decreasing L. For example, 8is increases 
from 52 at L = 4 through 56, 60 and to 70 as L decreases through 3, 2 and 1. 

The number of lattice knots of minimal length in Sl are displayed in tables |5] and 
[8] (for some compound knots) in the appendix. We display some of these results in 
figure |4l where we plot the number of minimal length lattice knots against the number 
of minimal length lattice knots in Si for different knot types on a log-log scale. The 
data scatter in the plot, showing that the number of minimal length lattice knots may 
change significantly with a decrease in L. 

Lattice knots are partitioned into symmetry classes due to invariance under 
rotations or (in the case of amphichircal knots) reflections which respect the orientation 
of Sl. Thus, the sets of lattice knots enumerated in tables |6] and |8] partition into 
symmetry classes. These classes are listed in tables |9] and [TO] in the appendix. The 
symmetry classes are denoted by 2 a 4 f, 6 c 8 d 12 e 16-' for each set of lattice knots. For 
example, for the knot type 3+#3f in table H] and for L = 1, the symmetry classes 
are 2 3 4 36 8 285 meaning that there are 3 classes with two members (each member is a 
lattice knot), 36 classes with 4 members and 285 with 8 members. 

The data show that the entropy decreases with decreasing L, if he the minimum 
length of the polygon does not change. In cases where the minimum length increases 
with decreasing L, it may be accompanied by a large decrease or increase in entropy. 
See for example the data for 8^" and 8g at L = 1 and L = 2 in table |H 

In what follows, we will use the data in these tables to determine the response 
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Figure 4. A scatter plot of the number of minimal length lattice knots p ni K (K) 
in Si (X-axis) and the number of minimal length lattice knots p nK (K) (Y-axis) 
on a log-log scale. 



of the lattice knots when forces are applied to squeeze them in slabs with hard walls, 
or in a model where the forces are mediated via a flexible membrane to the highest 
vertices in the lattice knots. 

3. A grafted lattice knot between hard walls 

The free energy of grafted lattice knots in a slab of width L is given by 

Tl = Energy - T ■ Entropy. (2) 

The lattice knots have fixed length (this is the canonical ensemble), and we assume in 
this model that the length is fixed at the minimal length in the slab Sl- 

In this model, the entropy should be given by logp% L K {K), where Ul,k is the 
minimal length of the lattice knot of type if in a slab of width L. 

We assign an energy to the lattice knot as follows: Compressing a minimal length 
lattice knot in a slab of width L will generally reduce its entropy, but at a minimum 
value of L, no further compression can take place because a minimal length lattice 
knot cannot be realised in a narrower slab. Instead, further increase of pressure on 
the Sl will induce forces along the edges of the polygon, and at a critical value of the 
force, these induced forces will overcome the elastic or tensile strength of the edges 
composing the lattice knot. The result is that the lattice knot will either stretch in 
length to fit in a narrower slab, or it will break apart and be destroyed. We assume 
the former case (our data show that in most cases the level of stretching is less than 
10% of the rest length of the lattice polygon). Thus, assign a Hooke energy to the 
polygon, with rest length equal to uk- The energy is then given by 

Energy = k ■ (ul,k - n K ) 2 , (3) 
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where til k IS the minimum length to accommodate the lattice knot in a slab of width 
L. There is no Hooke energy in the event that til,k — n K- 

With the above in mind, we define the free energy as follows 

Tl = k(n LiK - n K f - T log p^ L K (K). (4) 

For example, in the case of the unknot Oi, one may determine directly that — 4 and 
p noo = 1 while p ni = 3. This shows that J-q = — Tlog 1 = and T\ = — Tlog3. 

It is important to note that this free energy is a low temperature or a stiff Hookc 
spring approximation. It is low temperature because thermal fluctuations in the length 
of the lattice knot are not modeled (that is, the lattice knot is always in the shortest 
possible conformations in §^), and it is a stiff Hooke spring approximation for that 
same reason: the energy barrier to stretch the polygon to iij, ^ + 2 edges in length is 
too big, and those states do not make a measurable contribution to the free energy. 

Free energy differences as a result change in entropy and the Hooke term induces 
entropic forces pushing against the walls of the slab. These forces should push the 
walls apart, both decreasing the length of the lattice knot and increasing its entropy. 
They are given by 

F L = A X T L =T L - (5) 

If an externally applied force / squeezing the walls of together exceeds (that is, 
if |/| > then the force / will overcome the entropic and Hooke terms in the free 

energy and squeeze the lattice knot into Sl-i- Thus, the critical values of an applied 
force pushing against Fl are given when 

h = -F L (6) 

and if |/| > |/l| then the walls of Sl will be squeezed together to compress the lattice 
knot. We call /^ the critical force of the model. 

Compressing the lattice knot between two hard walls performs work on the knot, 
and conversely, if a lattice knot is placed in a narrow slab and the slab expands as a 
result, then the lattice knot performs work on the walls of the slab. The maximum 
amount of useful work that can be extracted from this process is given by 

L 

assuming that the expansion is isothermic and reversible. In our geometry, SL = 1. 
Thus, Wk reduces to J2l f L - 

For example, compressing a minimal length unknotted polygon between two hard 
walls reduces the entropy of the unknot polygon only when L transitions from 1 to 0. 
That is, the top wall in figure [1] can be pushed down without encountering resistance 
until L = 1. Further compressing to L = reduces the free energy by 

Fa = AxJx = {F 1 - Jo) = -Tlog3. (8) 

This shows that the critical force is /i = Tlog3 in this model. In this case, the unknot 
lattice polygon can do at most Wo 1 = T log 3 units of work, assuming that it is placed 
in §o arL d allowed to expand the slab. 
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3.1. Squeezing minimal length lattice trefoils between two planes 

These ideas can be extended to lattice knots, using the data in tables [5] and [5] 

In the case of the trefoil knot 3f, it follows from table [5] that n L 3 + = 24 if L > 2 
but that n 1 3 + = 26. In other words, the length of the lattice knot increases from 24 
to 26 if it is squeezed by a force into a slab §i . This stretching of the polygon stores 
work done by the compressing force in the form of elastic energy, which we indicate 
by the Hooke term in equation (j4|). 

If the compressing force is removed, the lattice knot will rebound to length 24, 
and expand the slab to width L = 2, performing work while doing so. In L = 2 
the lattice knot is not stretched, but it still suffers a reduction in entropy. Further 
expansion of the slab width to L = 3 increases entropy, and there is thus an entropic 
force pushing the hard walls apart until the lattice knot enters a state of maximum 
entropy for values of L large enough. 

With this in mind, the free energies of minimal length lattice knots of type 3^~ as 
a function of L may be obtained from the data in tables [5] and [51 The results are 

J"l(3^) = 4fc-Tlog36; 
Wf) = ~T log 152; 
J- 3 (3+) = -T log 1660; 
■F> 4 (3^) = -T log 1664. 

Observe the steady increase in the entropy of the lattice knot with increasing L. At 
L = 2 the knot cannot be compressed to L = 1 without increasing its length to 26 
edges, and the Hooke term Ak appears in .Fi(3^~). 

From the above data one may compute the critical forces for lattice knots of type 
3+. These are 



oo, if L = 1 

Ak + T log (38/9), ifL = 2 

T log (415/38), ifi = 3 

T log (416/415), ifi = 4 

0, if L > 5 



(9) 



The results for above show that there is no entropy loss with decreasing L until 
L = 4. Thereafter, compressing the lattice knot results in entropy loss, and the critical 
forces are given above. At L — 2 the knot stretches to accommodate conformations 
in L = 1, with the result that a Hooke term appears. Observe that fx — oo, since 
a non-trivial lattice knot cannot be squeezed into §0: unless the compressing force 
overcomes the strength of the edges and breaks the polygon apart. 

The above expressions for the critical forces gives a compression profile for 
squeezing the lattice trefoil between two planes. We illustrate this profile as a bargraph 
in figure [5] - where we put k = 1/4 and T = 1 . In this case one may also compute 
W 3 + = Ak + T log(416/9). The choice of k = 1/4 in figure [5] gives a Hooke contribution 
of one to the free energy if the polygon should stretch by two edges. At this level, the 
Hooke energy does not dominate the free energy. 

Similar data can be obtained for the figure eight knot 4i, and its critical forces 
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Figure 5. A profile of critical forces for the lattice knot of type St. 
Compressing the knot between two plates encounters no resistance for L > 5. 
If L = 4, then there is a small resistance (not visible on this scale), and for L < 3 
a larger resistance. Bars in red indicate that the critical forces are due to entropy 
reduction alone, and blue bars denote a Hooke contribution to the critical force. 
In this example, k = 1/4 and T = 1. 
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Figure 6. A profile of critical forces /j, for the lattice knot of type 4i. 
Compressing the knot between two plates encounters no resistance for L > 5. 
If L = 4, then there is a small resistance, and for L < 3 a larger resistance. 
Bars in red indicate that the critical forces are due to entropy reduction alone, 
and blue bars denote a Hooke contribution to the critical force. In this example, 
fci = k 2 = 1/4 and T = 1. 



are give by 



oo, if L = 1 

32fc + T log (758/185) , if L = 2 

4fc-Tlog (379/170), if £ = 3 

T log (114/85), ifL = 4 

0, if L > 5 



(10) 
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Figure 7. A profile of critical forces /j, for lattice knots of types 5j and 5j • 
There is no resistance to compression until L = 5. There are small resistances if 
L = 4, and this decreases even further for L = 3 in the case of 5]*~. In the case 
of 5^" the resistance increases with decreasing L. Bars in red indicate that the 
critical forces are due to entropy reduction alone, and blue bars denote a Hooke 
contribution to the critical force. In this example, k = 1/4 and T = 1. 



displayed in figure [6] 

In this knot type, the lattice knot increases in length both in the transition from 
L = 3 to L = 2, and then again to L = 1, as seen in table [5] While one will generally 
expect Jl > /l+i (that is, a larger force is necessary to compress the knot in narrower 
slabs), there is an interplay between entropy and the Hooke terms in the free energy, 
and in some cases /l < /l+i- This is for example seen in figure[6]in the data for L = 3 
and L = 4 for the knot type 4i. One may also verify that = 36fc + Tlog(456/185). 

The results for five crossing knots 5f and six crossing knots are illustrated in 
figures [7] and [8] 

Our numerical data on seven and eight crossings knots in tables [5] and [6] were used 
to compute critical forces and W for each of those knot types. Data on compound 
knots up to eight crossings are listed in tables [7| and [8j and we similarly determined 
critical forces and W for those knot types. The results are shown in tables Q] and [H 

3.2. Discussion 

Compressing a lattice knot between two hard walls decreases the entropy of the knot, 
until the walls are close enough together. Then the lattice knot expands laterally 
and in length as it finds conformations which can be accommodated in even narrower 
slabs. 

For example, a lattice trefoil loses entropy in Sl as L is reduced from L = 4 
to L = 2, but then has to increase in length by 2 if it is compressed into Si. This 
stretching of the lattice knot to a longer length changes its entropic properties, and 
may even increase its entropy, as the longer lattice polygon may be able to explore more 
states. However, the Hooke energy involved in stretching the lattice knot increases 
the free energy, and also increases the critical force necessary to compress the knot 
into a narrower slab. 
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Figure 8. A profile of critical forces fa for the lattice knots of types 6^, 
and 63. There is no resistance to compression until L = 5 in all cases. There 
are small resistances if L = 4. The negative bar between L = 2 and L = 3 for 
6]*~ shows that a gain in entropy overwhelms the Hooke forces when the lattice 
knot is compressed from L = 3toL = 2-in fact, no force is necessary, as the 
entropic force pulls the walls together. Bars in red indicate that the critical forces 
are due to entropy reduction alone, and blue bars denote a Hooke contribution to 
the critical force. In this example, k = 1/4 and T = 1. 
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Figure 9. A profile of critical forces fa for the compound lattice knots of types 
3^#3^ and 3]*~#3J~. There is no resistance to compression until L = 5 in all 
cases. Bars in red indicate that the critical forces are due to entropy reduction 
alone, and blue bars denote a Hooke contribution to the critical force. In this 
example, k = 1/4 and T = 1. 
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Figure 10. The same data as in figure [9] but with the vertical axes scaled 
logarithmically to enhance the data at larger values of L. 



The critical forces of the trefoil knot are listed in equation ©. These forces 
are induced by entropy loss when L > 2, but at L — 2 a Hooke energy contribution 
also appears. The sum of these forces gives the total amount of work in an isothermic 
compression of the lattice knot. We found this to be W 3 + = 4k + Tlog(416/9). 
There are two contributions to W 3 + , namely an entropic and a Hooke contribution. 
These contributions are equal in magnitude at a critical value of k (or equivalently, 
a critical value of T). In the case of the trefoil knot this critical value of k is 
fc 3+ = [Tlog(416/9)]/4 « (0.95836 . . .)T. If k > k c , then the Hooke term dominates 
W 3 + and if k < k c , then the work done has a larger contribution from entropy reduction 
in the process. 

In the case of the figure eight knot, one may similarly determine the critical value 
of k: k4 1 — (0.02505 . . .)T. In this case, k c is very small, and the total amount of 
work done in compressing the knot to L — 2 is dominated by the Hooke term even 
for relatively small values of k. The two five crossing knots have k 5 + — (0.23974 . . .)T 
and k 5 + = (0.06167 . . .)T. 

Data on other knots are listed in tables [1] and [2] The critical values of k can be 
determined by solving for k from the last column in these tables. Observe that some 
knot types have negative values of k c , for example, k e + = — 0.015025T. This implies 
that a large entropy gain occurs when the knot stretches in length to fit in Si, and 
this can only be matched by a negative Hooke constant. 
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Table 1. Critical forces Jx for knots up to seven crossings and for compound knot types up to eight crossings. The last column is the 
maximum amount of useful work which can be extracted if the lattice knot expands against the hard walls of the slab, pushing them apart 
from L = 1. In all these cases we set T = 1. 
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Table 2. Critical forces for knots with eight crossings. The last column is the maximum amount of useful work which can be extracted 
if the lattice knot expands against the hard walls of the slab, pushing them apart from L = 1. In all these cases we set T = 1. 



Squeezed Lattice Knots 



16 



4. A grafted lattice knot pushed by a force 

In this section we consider the model inspired by figure [2j An external force / 
compresses a grafted lattice knot onto a hard wall. If p n (h,K) is the number of 
such lattice knots with highest vertices at height h above the bottom wall and length 
nh-K, then the partition function is given by 

oo 

Z(f) = £ PnhiK (h,K)e- E *-f h (11) 

n,h=0 

where rih-x = rnin^<^{nL : x} is the minimal length of the polygons in §l for all L < h. 
For example, %. 3 + = 26 and n h . 3 + = 24 for all h > 1. 

In the partition function, positive values of / mean that the vertices in the top 
layer are pushed towards the bottom wall, and negative values / mean that the force 
is pulling the vertices in the top layer from the bottom wall. Observe that we put the 
Boltzman factor ksT = 1 in this definition, and so use lattice units throughout. 

The function Eh is an energy of the lattice knot. We shall again use a Hookc 
energy for the polygons, namely 

E h = k{n h ., K - n K f . (12) 

In addition, we assume a low temperature approximation, namely that only polygons 
of the shortest length in §l contribute to Z(f) for all L. That is, in the case of 3^~, 
only polygons of lengths 24 and 26 (when h = 1) contribute. This approximation is 
also valid in the regime that the Hooke constant k is large. 

In other word, the low temperature and large Hooke constant approximation of 
the partition function is 

oo 

Z*(f) = Y,Pn h , K (h,K)e- k ^- n ^ 2 -f h . (13) 
h=a 

For / < we observe that the polygon is pulled by its highest vertices from the bottom 
wall, and that it will stretch in length if the forces overcome the tensile strength of the 
edges. Stretching the edges in this way will cause the Hooke energy term to increase 
quadratically in n oc h, while the force / couples only linearly in with h. Thus, this 
regime will be a purely Hooke regime, provided that k is large enough. 

Thus, we obtain a model of fixed length lattice knots, which may stretch to longer 
states if pushed against the bottom wall by large forces to accommodate itself into a 
conformation with small h. 

The (extensive) free energy in this model is given by 

J)=logZ*(/) (14) 

and its derivatives give the thermodynamic observables of the model. For example, 
the mean height of the grafted lattice knot is 

_ dT ! j:Z hp n (K;h)e- k ^-^ 2 -f h 

(ft,K "T" ' (15) 

while the second derivative 

df (h) K df 

is the fractional rate of change in mean height with /, and is a measure of the linear 
compressibility of the lattice knot due to a force acting on its highest vertices. 
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Figure 11. The mean height {h) (red curve) and compressibility ko 1 (blue 
curve) as a function of / for the unknot 0i. For negative values of / (pulling 
forces) the mean height is 1, and as the lattice knot is compressed into pushing 
forces, its mean height decreases until it approaches zero. Note that (h) = 2/3 if 
/ = 0. 



The data in tables [5] and |6] can be used to determine kk- For example, for the 
unknot grafted to the bottom wall, one has n^-K = nL,K — 4 if L = h = or 
L = h — 1. Thus, we determine the partition function in this case to be 

Z* i (f) = l + 2e-f (17) 

since Eh = for both h = and h = 1 in this case. Observe that P4(0, 0i) = 1 and 
P4(l, 0i) = 2 in this model. 

One may now compute the mean height and kq 1 for the unknot directly from the 
above. The results are 

w * = iHf7' K0i = irl^7 (18) 

and these are plotted in figure 111! If this lattice knot is released from its maximum 
compressed state in a slab of width L = 0, and allowed to expand, then work can be 
extracted from the expansion. The maximum amount of work which may be extracted 
is given by the free energy differences between the L = slab and the free energy at 
/ = 0. This is given by 

Woi = T } | f=Q = T f | L=0 = log 3 - log 1 - log 3. (19) 

This is the same value obtained in the first model. 

4--1- Compressing minimal length lattice trefoil knots 

By noting that n h . 3 + = 24 if h > 2 and ^.3+ = 26 in table O one may determine 
Pn(h 3+)(^> ^1") by examining the data in table [SJ 

In particular, it is apparent that p n M 3 +j (1, 3f) = 36 and p n , 2 3 +)(2, 3^ ) = 152. 
However, these 152 lattice knots of length n = 24 in S2 are also counted in §3, and 
so must be subtracted from the data in column L = 3 to obtain p n ^ 3 3 +-j(3,3^). In 

particular, it follows that p n ^ 3 3 +)(3, 3j ) = 1660 — 152 = 1508. 

Similarly, one may show that p n ( 4 . 3 +) (4, 3^ ) = 4 and p n ^ >5 . 3 +~ ) (> 5, 3^) = 0. This 
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Figure 12. The mean height (h) + and compressibility « + of grafted lattice 

knots of type 3^. Negative forces are pulling forces, stretching the lattice knot 
from the bottom plane. Positive forces are pushing forces. The mean height 
decreases in steps from a maximum of about 4 to about 3 at / = and then to 
1 for large positive /. There are two peaks in ft + . The highest peak at positive 

f corresponds to the pushing force overcoming the Hooke term in J 7 , increasing 
the length of the lattice from 24 to 26 and pushing it into a slab of width L = 1. 
In this example, k = 1/4. 



shows that 

Z*(3+) = 36e- 4k - f + 152e- 2f + 1508 e~ 3/ + 4 e~ 4/ . (20) 
The mean height of the lattice knot is 

9e- 4fc + 76e- / + 1131e- 2/ +4e- 3 - f 

( ' 3 t ~ 9 e~ 4 k + 38 e- f + 377 e~ 2 f + e" 3 f { ' 

The (linear) compressibility of the lattice knot of type 3^~ is a more complicated 
expression, given by 

(342 + 13572 e~f + 81 e~ 2/ ) er ik -f + (14326 + 152 e^ + 377e" 2/ ) e^ 3/ 
Ks t = (9e- 4fc + 76e~f + 1131 e~ 2 / + 4e" 3 /) (9e" 4fc + 38 e~f + 377e~ 2 f + e~ 3 f)' 
In figures [T2l and [TBI the mean height and k 3 + are plotted as a function of the force for 
k = 1/4 (figure [T2|) and k — 4/3 (figure fT3|). In figure [12] the mean height decreases 
to L = 1 in two steps, the first at negative (pulling) forces, and the second a step 
from a height of roughly h = 3 to h = 1 . In this step (which shows up as a peak in 



figure 12(b) ) the polygon is squeezed into Si as the force overcomes both the entropy 
reduction and Hooke term. 

Increasing the Hooke constant k produces graphs similar to figure [131 There 
are now three peaks in k, each corresponding to a reduction of the lattice knot from 
height h to height h — 1 as the applied force first overcomes entropy and then the 
Hooke energy to push the lattice knot into Si . 

The total amount of work done by letting the lattice polygon expand at zero force 
from its maximal compressed state in L — 1 is given by W 3 + = J~j\ L=1 — Tf\^_Q. 
This gives 

W 3 + = log (1664 + 36e" 4fc ) - log (36e- 4fe ) = log (*|S e 4fc + l) . (22) 
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Figure 13. The similar plots to figure Il2l but with k = 4/3. This larger Hookc 
energy requires a larger force to push the lattice knots from a L = 2 slab into Si- 
This shows up as a third peak in the ft-graph above. 



4-2. Compressing minimal length lattice figure eight knots 

The minimal length of a lattice polygon of knot type 4i (the figure eight knot) is 36 
in Si, 32 in §2 and 30 in Sl with L > 3. In other words, states with heights 1 or 2 
will have a Hooke energy, as they have been stretched in length to squeeze into slabs 
with small height. 

By consulting the data in tables [5] and [6l the partition function of this model can 
be determined, and it is given by 

Z*(4i) = 1480e- / - 36fe + 6064e- 2/ - 4fe + 2720e- 3/ + 928e~ 4/ . (23) 

The mean height of this lattice knot is given by 

_ I85e- 36k + 1516e-- / '- 4 ' £ + 1020e- 2/ +464e" 3 ^ 
' ' 4l ~~ 185e~ 36A; + 758e-/- 4fc +340e~ 2 / + 116e" 3 -f 



and it is plotted in figure 14(a) for k — 1/4. Observe that the mean height decreases 
in steps with increasing /, from height 4 to height 2, before it is squeezed into Si for 
sufficiently large values of /. 

The compressibility is plotted in figure [l4(b)| and the two peaks correspond to the 
decreases in (/i) 4 in steps. At the peaks, the lattice knot has a maximum response 
to changes in /. The expression for the compressibility is lengthy and will not be 
reproduced here. 

Finally, the total amount of work done by releasing the knot from Si at zero 
pressure and letting it expand isothermically, is given by 

W4l =log(f|e 36fc + |ie 32fc + l) (25) 

In the case that k = 0, W 3 + « 3.855 > 2.023 re W^. In other words, more work 
is done by the trefoil knot. However, if k = 1/4 (and a Hooke term is present), then 
the relation is the opposite: W 3 + re 4.841 < 10/379 re W^. Equality is obtained 
when k = 0.06588 .... 
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Figure 14. The mean height (/i) 4l and compressibility K4 1 of grafted lattice 
knots of type 4i . Negative forces are pulling forces, stretching the lattice knot 
from the bottom plane. Positive forces are pushing forces. The mean height 
decreases in steps from a maximum of 4 to about 2.8 at / = and then to 1 
for large positive /. There are two peaks in K4 ± . Both peaks correspond to a 
reduction in the slab width L as the force first overcomes the Hooke term from 
L = 3 to L = 2, and then again at larger pushing forces, the Hooke term from 
L = 2 to L = 1. In this example, k = 1/4. 



4-3. Compressing minimal length lattice knots of types 5j and 5^ 

The partition functions of minimal lattice knots of types and are given by 

Z*(5f) = 72e' f ' 16k + 760e _2/ + 600e~ 3/ + 1936e _4/ +40e~ 5/ ; 

Z*(5j) = 6240e~ / ~ 36fc + 720e~ 2/ + 24072 e~ 3/ + 30544 e~ 4/ + 2120 e~ 5 f . (26) 

These expressions show that the maximum height is 5 while there are Hooke terms 
for the transition from h = 2 to h = 1. 

The mean heights can be computed, and are given by 
9e -i6fe + i9 0e -/ + 225e~ 2 ^ + gOSe-^' + 25e" 4/ 



(h) 



5 t 9 e- 16 k + 95 e~f + 75 e- 2 / + 242 e" 3 f + 5 e" 4 / 

_ 780e~ 36fc + 180 e-f + 9027 e~ 2f + 15272e~ 3/ + 1325e~ 4/ 

^ ' 5 i ~ 780e" 36fc + 90 e-/ + 3009 e-2/ + 3818 e- 3 / +265 e- 4 / ( ' 



Wk are similar given by 



W 5 + = log (if e ib,£ + 1) and W s + = log e 3bfc + l) . (28) 

The results for the 5-crossing knots are plotted in figure [TS] The curves for the 
mean height (h) are very similar , with 5^~ undergoing a smoother compression with 
increasing /. The knot type 5^ shows more variation, and this is very visible in 



the plots for k in figure 15(b) where there are several sharp peaks for 5^, but less 
pronounced changes for 5 J. In both knot types, the Hooke constant is k = 1/4. 

4-4- Compressing minimal length lattice knots of types 6^, and 63 

The results for 6-crossing knot types are displayed in figure [16] These knot types 
exhibit more similar behaviour than the two 5-crossing knot types. 
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Figure 15. The mean height (/i) 5 and compressibility Ks„ of grafted lattice 
knots of types 5^" (red curves) and 5 J (blue curves). The mean height decreases 
in steps from a maximum of 5 to 1 for large positive /. The transition for is 
smoother than for 5^, which exhibits peaks in k corresponding to critical forces 
overcoming the Hooke terms and entropy. There is one low and are two prominent 
peaks in n +. In this example, k = 1/4. 





30 



/ / 

(a) The mean height of 6-crossing knots (b) The compressibility of 6-crossing knots 



Figure 16. The mean height {h)g and compressibility «6* of grafted lattice 
knots of types 6^" (red curves), 6 J (blue curves) and 63 (green curves). The mean 
height decreases in steps from a maximum of 5 to 1 for large positive /. These 
knot types follow a very similar pattern, with 63 posing the most resistance to 
compression. Each knot also exhibits two peaks in k of similar height. In this 
example, k = 1/4. 
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Table 3. The work done by releasing the lattice knot in §1 and letting it expand 
to equilibrium at zero pressure. 



4-5. Discussion 



The compression of lattice knots typically show a few peaks in the compressibility. 
Since the compressibility is the defined as the fractional change in the width of the 
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Table 4. The work done by releasing the lattice knot in §i and letting it expand 
to equilibrium at zero pressure. 

lattice knot with incrementing force, a high compressibility corresponds to a "soft" 
lattice knot, and a low compressibility to a "hard" lattice knot. The peaks observed 
in figures from figure [T2] to figure [TB] correspond to values of / where the lattice knots 
are soft. 

Our general observation is that the compressibility of lattice knots in this model is 
not monotonic: That is, the lattice knot does not become increasingly more resistant 
to further compression with increasing force. Instead, it may alternately become 
more or less resistant to compression, as the compressibility goes through peaks ( "soft 
regimes" ) and troughs ( "hard regimes" ) . 

The work done by releasing the lattice knot in §i and letting it expand at zero 
pressure was computed above for knots up to 5-crossings. For example, if the lattice 
unknot is placed in §o an d allowed to expand isothermically to equilibrium at / = 0, 
then Wqj = log 3. 

A similar calculation gave W 3 + = log (^pe 4 * + l) for expansion from Si, and 
expressions were also determined for other knot types. The results are given in tables 
[3] and |4] for the remaining knot types. 

Observe that these expressions are different from those in the model in section 3: 
There the data were obtained by computing the sum over all the critical forces. In each 
case the minimal values of the forces required to overcome resistance by the squeezed 
polygon were determined. This is a different process to the above, where the lattice 
knots are placed in Si and then allowed to expand isothermically to equilibrium. The 
total work was obtained by computing the free energy difference between the initial 
and final states. 

Finally, we examined the effect of the Hooke parameter k by comparing W for 
two different values of fc, namely fc = 1/4 and k = 1. For each knot type, the point 
(Wfc = i/ 4 , Wfe=i) is plotted as a scatter plot in figure [T71 These points line up along 
a band, indicating that the Hooke term makes a significant contribution to the free 
energy. 

5. Conclusions 

In this paper we presented data on minimal length lattice polygons grafted to the 
bottom walls of slabs Sl- We collected data on these polygons by implementing the 
GAS-algorithm for knotted lattice polygons in a slab in the lattice, and sieving out 
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Figure 17. A scatter plot of W with fc = 1/4 (horizontal axis) and k = 1, (vertical 
axis). The data for prime knot types accumulate along a band, indicating that 
the Hooke term makes a significant contribution to the free energy. 



minimal length polygons. Our data were presented in tables [5] and [6] for prime knot 
types, and tables [7] and [8] for compound knot types. Our results verify, up to single 
exceptions, the data obtained in reference [8]. 

We determined the compressibility properties of minimal lattice polygons in two 
ensembles. The first was a model of a ring polymer grafted to the bottom wall of a 
slab with hard walls, and squeezed by the top wall into a narrower slab. 

The second was a model of a ring polymer grafted to a hard wall, and then pushed 
by a force / towards the bottom wall. 

In these models we determined free energies, from which thermodynamic 
quantities were obtained. In the first model we computed the critical forces which 
squeeze the lattice polygons into narrower slabs, and we computed the total amount 
of work that would be done if the lattice polygon was squeezed into Si . These data 
are displayed in tables Q] and [2j 

The profiles of critical forces in this model are dependent on knot types, and 
generally increases with decreasing values of L. The dependency on knot type shows 
that entanglements in the lattice knot, in addition to the Hooke term and the entropy, 
plays a role in determining the resistance of the lattice knot to be squeezed in ever 
thinner slabs. 

We examined the mean height and compressibility profiles for some simple knot 
types in section 4. The amount of work done by expanding expanding lattice knots 
from maximal compression at zero force was also determined, and the results were 
listed in table [3] and |4j 

In contrast with the results in section 3, we computed the compressibility of 
lattice knots as a function of the applied force in section 4. Our results show that 
the compressibility is dependent on the knot type and the Hooke energy and is not 
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a monotonic function of the applied force. In many cases there are peaks in the 
compressibility where the lattice knot becomes (relatively) softer (more compressible) 
with increasing force. This may feature may be dependent on either the entanglements, 
or the loss of entropy, or the Hooke term, or on a combination of these. These 
observations indicate that such peaks may be seen in some conditions when knotted 
ring polymers are compressed. 

Finally, the inclusion of other energy terms, for example a bending term, or a 
binding term to to the walls of the slab, can be done. We expect such additional 
terms to have an effect on our results. For example, a bending energy will make the 
lattice knot more rigid and less compressible, increasing the critical forces computed in 
section 3, and reducing the heights of the peaks seen in the compressibility in section 4. 
A further generalisation would be to simulate lattice knots in other confined spaces, 
for example in pores or in channels. Such models pose unique questions, since the 
algorithm is not known to be irreducible. 
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Appendix: Numerical Results 

Our raw data on the estimates of the exact minimal lengths and number of minimal 
length lattice knots in §l are displayed in tables ??, [71 El and[5] 

Data on the symmetry classes of lattice knots confined to slabs §£. Data are 
presented in the format A a B b C c . . . denoting a symmetry classes of polygons, each 
with A members, b symmetry classes of polygons, each with B members, and so on. 
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Table 5. The minimal length of lattice knots of prime knot type to eight crossings 
confined to slabs of width L. Cases with an increase in minimal length are 
underlined. 
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Table 6. The number of lattice knots of prime knot type to eight crossings 
confined to slabs of width L. 
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Table 7. The minimal length of lattice knots of compound knot type to eight 
crossings confined to slabs of width L. Cases with an increase in minimal length 
are underlined. 
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Table 8. The number of lattice knots of compound knot type to eight crossings 
confined to slabs of width L. 
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Table 9. The symmetry classes of prime lattice knots upto seven crossings, and compound lattice knots upto eight crossings, confined to 
S L with L < 8. 
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24 1910 


24 1910 


24 1910 


83 


4 3gll35 


g250 


8 1 


8 1 


12 1 


12 1 


12 1 


12 1 


8 4 + 


4 63g2493 


gll90 


g230 


gl47 16 831 24 19 


16 58 24 939 


24 997 


24 997 


24 997 


si 


2 2 4 5 8 6 


^57g50994 


4 10 g 613 


16 24 


16 24 


16 24 


16 24 


16 24 


8 6 + 


gl205 


8 1 


g3013 16 5 


16 230 


24 230 


24 230 


24 230 


24 230 


8+ 


g2429 


g6172 


8 1 


16 1 


24 1 


24 1 


24 1 


24 1 


8 8 + 


gl696 


g347 


g4959-^g38 


16 65 


24 65 


24 65 


24 65 


24 65 


89 


gl908 


g2130 


g 239 16 134 24 36 


gl 16 1488 24 282 


12 1 24 1770 


12 1 24 1770 


12 1 24 1770 


12 1 24 1770 


8 io 


g26 


g5964 


g663 


16 35 


24 35 


24 35 


24 35 


24 35 


si 


g 488 


g453 


g5056 


16 4 


24 4 


24 4 


24 4 


24 4 


812 


gl830 


g2S970 


g36 


16 108 


24 108 


24 108 


24 108 


24 108 


813 


gl9974 


g!90 


glO 


gl21 16 316 24 107 


16 28 24 S16 


24 544 


24 544 


24 544 


8 14 


g7375 


g860 


g748 16 ll 


16 15 


24 15 


24 15 


24 15 


24 15 


815 


4 2 8 x 


gllO 


4 10g559 


4 7g350 12 3 16 1259 24 57 


12 10 24 1666 


12 10 24 1666 


12 10 24 1666 


12 10 24 1666 


8 16 


4 14g4083 


g!46 


8 6 


16 2 


24 2 


24 2 


24 2 


24 2 


8l7 


g7542 


gl036 


g3234 


gll54 16 998 24 64 


16 60 24 2156 


24 2216 


24 2216 


24 2216 


8l8 


4 4 8 2 


2 6 4 6g90 


g8500 


g!48 


24 148 


24 148 


24 148 


24 148 


819 


1 l 2 3 4 9gl5 


g207 


4 13g75 16 30 


g 12 12 13 16 30 24 249 


12 25 24 279 


12 25 24 279 


12 25 24 279 


12 25 24 279 


8 20 


4 3 8 13 


g!59 


8 5 


16 6 


24 6 


24 5 


24 5 


24 5 


8 2 + i 


4 2g2 


g90 


4 lg428 


12 1 16 315 24 852 


12 1 24 1167 


12 1 24 H67 


12 1 24 H67 


12 1 24 1167 



to 



Table 10. The symmetry classes of prime lattice knots up to eight crossings, confined to Sl with L < 8. 



